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Abstract 

We report preliminary results on the measurement of branching fractions of exclusive decays of 
neutral and charged B mesons into two-body final states containing a charmonium state and a 
light strange meson. The charmonium mesons considered are J/V', i^C^S) and Xci, and the light 
mesons are either K or K*. We use a sample of about 124 million BB events collected with the 
BABAR detector at the PEP-II storage ring at the Stanford Linear Accelerator Center. 
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1 Introduction 



Fully hadronic decays of B mesons have proven to be an effective laboratory to study and provide 
tests of the theory of heavy quarks as well as the dynamics of strong interactions in heavy meson 
systems. The tree level diagram of the decays under study is shown in Figure ^ 




Figure 1: Tree level diagram of a meson decaying into a charmonium state and a kaon. 

The dynamics of the decay is expected to be highly affected by strong interactions effects, 
especially by the long distance non-pertubative aspect of QCD. There are various phenomenological 
approaches to treat these decays, which provide different estimates for the branching fractions (see 
^ and references [2-12] therein). 

Charge asymmetry measurements can be a powerful tool for seeking new physics. The Standard 
Model predicts small direct CP violation 2 , thus large charge asymmetries would indicate new 
physics 0. 

The list of the branching fractions measured and decay modes considered in this paper is shown 
in Tabled 

Table 1: Branching fractions and decay modes considered in this analysis. 



Decay Channel 


Secondary decay mode 
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J/iP K*^ 


K*^' - 








m - 




B+ - 


J/ip K*+ 


K*+ - 












B^ - 












vrO ^ 
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B+ - 


J 111) K+ 














B^- 


ip{2S) K*^ 


K*^ - 


K+ 






ij{2S) 


^ £+£- 


B+ - 


i){2S) K*+ 


K*+ - 




-A 








B^ - 


■> i}{2S) Kl 
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B+ - 


il){2S) K+ 














B^ - 


-> Xci K*^ 


K*^ - 


K+ 


TT", 




Xci 




B+ - 


- Xci K*+ 


K*+ - 




-A 




J/^ - 


e+£- 


5° - 


- Xci 














B+ - 


- Xci K+ 












77 
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2 The BaBAR detector and dataset 



The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e'^e~ storage ring from 1999 to 2003. This represents a total integrated luminosity of 112.4 fb~^ 
taken on the T(4S) resonance, producing a sample of 123.95 million BB events. 

The BABAR detector is described elsewhere jl]. Surrounding the interaction point, a 5 layer 
double-sided silicon vertex tracker (SVT) provides precise reconstruction of track angles and B 
decay vertices. A 40 layer drift chamber (DCH) provides measurements of the transverse momenta 
of charged particles. An internally reflecting ring-imaging Cherenkov detector (DIRC) is used 
for particle identification. A CsI(Tl) crystal electromagnetic calorimeter (EMC) is used to detect 
photons and electrons. The calorimeter is surrounded by a 1.5T magnetic field. The flux return is 
instrumented with resistive plate chambers (IFR) used for muon and neutral hadron identification. 

3 Analysis Method 

Multihadron events are selected by demanding a minimum of three reconstructed charged tracks 
in the polar-angle range 0.41 < Oiab < 2.54 rad. Charged tracks must be reconstructed in the DCH 
and are required to originate at the nominal beamspot, within 1.5 cm in the plane transverse to 
the beam and 10 cm along the beam. Events are required to have a primary vertex within 0.5 cm 
of the average position of the interaction point in the plane transverse to the beamline, and within 
6 cm longitudinally. 

Charged tracks used in this analysis are required to include at least 12 DCH hits, to have a 
transverse momentum >100 MeV/c. 

Photons are reconstructed from EMC clusters. The radial energy profile (LAT) [3] of the 
cluster is used to discriminate electromagnetic from hadronic clusters. Photons are required to 
have a minimum energy of 30 MeV, a radial energy profile less than 0.8, and to be in the fiducial 
volume 0.41 < 9 < 2.41 rad. 

Electron candidates are selected using information from the EMC (radial energy profile and 
Zernike moment ^442 [HI ) , the ratio of the energy measured in the EMC to the momentum measured 
by the tracking system (E/p), energy loss (dE/dx) in the drift chamber and the Cherenkov angle 
measured in the DIRC. Electrons are also required to be in the fiducial volume 0.41 < 6 < 2.41 
rad. 

Muon candidates are selected using information from the EMC (energy deposition consistent 
with a minimum ionizing particle) and the distribution of hits in the IFR. Muons are required to 
be in the fiducial volume 0.3 < 6 < 2.7 rad. 

The charged kaon and pion candidates are selected using information from the energy loss in 
the SVT and DCH, and the Cherenkov angle measured in the DIRC. Kaon candidates are required 
to be in the fiducial volume 0.45 < 6 < 2.45 rad. 

The next step in the analysis is to combine tracks and/or neutral clusters to form candidates. If 
a particle decays through an intermediate state, this is constrained to its known mass, except for 
the K*. The selection has been optimized by maximizing the ratio S/\/S + B, where S and B are 
respectively the number of expected signal and background events obtained from GEANT4-based 
Monte Carlo simulation after the selection. 

The J/ip candidates are required to have an invariant mass 2.95 < Mg+g- < 3.14 GeV/cP and 
3.06 < M^+^- < 3.14 for J/jp e^e~ and J/ijj Ai"*"/^" decays respectively. 
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The '4>{2S) candidates are required to have an invariant mass 3.44 < Mg+g- < 3.74 GeV/c^ and 
3.64 < M^+^- < 3.74 GeV/c^ for 7p{2S) e+e" and ip{2S) ^ decays respectively. 

For J/ip — > e"^e~ and ip{2S) — > e'^e~ decays, electron candidates are combined with photon 
candidates in order to recover some of the energy lost through bremsstrahlung. Photons are required 
to be within 35mrad in polar angle from the electron track, and to have an azimuthal angle 
intermediate between the initial track direction (estimated by subtracting 50 mrad opposite to the 
bend direction of the reconstructed track) and the centroid of the EMC cluster arising from the 
track. 

In the Xci reconstruction {xd — > J/i^ 7), J/ip candidates are selected as described above. The 
associated 7 has to fulfill the following requirements: radial energy profile less than 0.8, Zernike 
moment ^442 less than 0.15 and energy greater than 0.15 GeV. Furthermore, Xci candidates are 
required to satisfy 0.35 < Mi+i-^ - Mi+i- < 0.45 GeV/c^. 

The vr'^ — > 77 candidates are required to satisfy 0.113 < M^^ < 0.153 GeV/c^. The radial 
energy profile of both photons are required to be less than 0.8. The energy of the soft photon has 
to be greater than 0.050 GeV and the energy of the hard photon has to be greater than 0.150 GeV. 

The 7r+ tt" candidates are required to satisfy 0.489 < M^+^- < 0.507 GeV/c^. The fol- 

lowing selection is also required: the Kg vertex has be more than 1mm from the charmonium vertex, 
and the angle in the x-y plane between the Kg momentum and the line joining the charmonium 
and Kg vertices has to be smaller than 0.2 rad. 

The K*^ and K*~^ candidates are respectively required to satisfy 0.796 < Mktt < 0.996 GeV/c^ 
and 0.792 < Mkw < 0.992 GeV/c^. In addition, for channels having a vr'' in the final state, the 
cosine of the angle between the K momentum defined in the K* rest frame and the K* momentum 
defined in the B rest frame has to be smaller than 0.8 (this helps in removing background coming 
from events with soft pions). 

Finally, B candidates are reconstructed by combining charmonium and kaon meson candidates 
and are selected by the use of two kinematic variables: the difference between the reconstructed 
energy of the B candidate and the beam energy in the center-of-mass frame AE = — E^^^^, and 

the beam energy substituted mass mES) defined as rriES = \J ^beam ~ P*B (^he * refers to quantities 
in the center of mass). For a true B meson, AE is expected to peak at zero, and the energy 
substituted mass ttt-es should peak at the B meson mass 5.279 GeV/c^. Only one reconstructed 
B meson is allowed per event. For events that have multiple candidates, the candidate having the 
smallest AE is chosen. Depending on the channel, around 10% of the candidates are removed by 
requesting a single B meson per event. The analysis is performed in the mES vs AE plane, defined 
as: 5.2 < mEs < 5.3 GeV/c^ and -0.12 < AE < 0.12 GeV. As an example, Figure [21 shows the 
AE and ttt-es distributions for the B J/ip K*^ (K~^ vr") channel. We subsequently define a 
signal box region in the rriEs vs AE plane, where the sensitivity is optimal. The signal box region 
is channel-dependent. For most of the channels, the signal regions are taken as the mean value 
±3(7 for both AE and m-ES- For channels with less statistics {'ip{2S) K* and Xci K* channels), the 
mES signal region was taken as 5.27 < uies < 5.29 GeV/c^, and the AE' signal region was taken 
as |A£^| < 0.04 GeV for channels with a vr'' in the final state and \AE\ < 0.03 GeV for the other 
channels. 

The selection efficiencies for each mode are obtained from Monte Carlo and are given by the 
number of expected signal events divided by the total number of generated events. While the 
Monte Carlo has been tuned to be as realistic as possible, one still has to correct for residual 
differences between data and simulated events. We have therefore applied additional corrections 
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Figure 2: AE and mss distributions for the B J/iJj K*^ {K^ vr^) channel. The blue points 
represent the Data and the histogram represents the Monte Carlo. An offset between the Monte 
Carlo and data distributions can be seen. It has been corrected. 



to the selection efficiency coming from particle identification, neutral particle, tracking, and Kg 
corrections. 

The number of signal events Ns is determined from the number of candidate events A'cand after 
subtracting the background. The mss distribution within the /S.E signal region is fitted by an 
Argus function [Jj and a Gaussian, and both functions are subsequently integrated within the 
mES signal region. The number of candidate events is given by the Gaussian integral. There are 
two components to the background: the combinatorial background and a peaking component (the 
component of the background that has /S.E and ttt-es distributions peaking at /S.E =0 GeV and m-Es 
=5.279 GeV/c^ respectively). The combinatorial background is obtained by integrating the Argus 
function within the mES and /S.E signal regions. The peaking component is obtained from Monte 
Carlo. There are two contributions to the peaking background. The first contribution is coming 
from feed-across events which, in the case of the J/ifj K*^ {K^ -k^) reconstruction, for instance, 
come from J/V' K*^ {K+ vr"), J/^jj K*+ {K^ and J/il) K*+ {K+ 7r°). The second contribution 
is coming from inclusive charmonium. For each of the contributions, the tties distribution is fitted 
within the IS.E signal region by an Argus and a Gaussian function, which are subsequently integrated 
within the tit-es signal region. The amount of peaking background is given by the Gaussian integral. 

The branching fractions are obtained from: 



BE = ^ (1) 

where is the number of BB events, e is the selection efficiency and / is the total secondary 
branching fraction. For channels with a K* in the final state, the feed-across contribution, which 
depends on the branching fractions that are being measured, to the peaking background can be 
important. Therefore an iterative procedure has been employed in which the feed-across contri- 
bution is re-estimated at each iteration. The procedure converges quickly as the feed-across is a 
small fraction of the number of signal events. When allowed by the size of the data sample, the 
branching fractions have been measured for both J/ij: — > e^e~ and J/^ — > decays separately. 
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Figure 3: mss distributions and fits within the AE Signal Box region for B charmonium K* 
channels. The top row represents the distributions for the J/ip K* channels, the middle row the 
ip{2S) K* channels, and the bottom row the Xci K* channel. From left to right, the columns show 
the distributions for the K*° K^tt^, K*^ K+ tt- , K*+ 7r+, and K*+ K+ vr" decay 

modes. The dashed lines show the combinatorial contribution to the background. The dotted lines 
show the peaking background contribution. 

The mES distributions within the AE signal region for candidate events are shown on Figures 
ElandU 

4 Systematic studies 

The systematic errors arise from the uncertainty on the number of BB events, the secondary 
branching fraction, the estimate of the selection efficiency, and the knowledge of the background. 

The systematic uncertainty on the number of BB events is 1.1%. It is common to all the 
branching fraction measurements. The secondary branching fractions and their errors have been 
taken from [Hj. 

For the tracking efficiency, we have used a flat correction of 0.8% per track with an associated 
error of 1.3% per track. The Kg efficiency corrections have been determined by the use of control 
samples and its errors from varying the Kg selection. The resulting error on the K^ efficiency 
varies from 0.8% to 2.0% depending on the channel. The uncertainty on the detection and energy 
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Figure 4: rnES distributions and fits within the AE Signal Box region for B charnionium Kg 
channels (top row) and B charmonium channels (bottom). From left to right, the columns 
show the distributions for the J/'ip,jp{2S) and Xci channels. 



Table 2: Breakdown of contributions to the systematic errors for the J/ip K* channels. All values 
are expressed relative to the measured branching fractions, in percent. 







K+ TT- 










(-'-(- irir 






B counting 


1.1 


1.1 


1.1 


1.1 


Tracking 


2.6 


5.2 


3.9 


3.9 


Polarizarion 


7.37 


3.85 


4.85 


6.86 




0.7 




0.8 




Neutral 


5.72 5.27 






5.14 5.61 


Second BF 


1.69 1.71 


1.69 1.70 


1.70 1.71 


1.69 1.70 


PID 


0.95 2.14 


0.88 3.57 


0.90 2.93 


0.61 3.73 


Background 


3.57 3.25 


1.18 0.91 


1.45 0.96 


2.34 2.04 


MC statistics 


1.24 1.23 


1.37 1.32 


1.51 1.45 


1.80 1.71 


Total 


10.65 10.49 


7.06 7.82 


6.97 7.42 


10.08 10.90 



measurement of photons is 2.5% common to all channels plus an additional channel-dependent 
correction. The uncertainty on the 7r° reconstruction is 5.0% for all channels plus an additional 
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Table 3: Breakdown of contributions to the systematic errors for the tp{2S) K* channels. All 
values are expressed relative to the measured branching fractions, in percent. 







K+ ir- 


Kl 7T- 




B counting 


1.1 


1.1 


1.1 


1.1 


Tracking 


2.6 


5.2 


3.9 


3.9 




2.0 




1.9 




Neutral 


7.2 






6.2 


Second BF 


11.74 


11.72 


11.72 


11.72 


PID 


0.97 


1.62 


0.21 


0.52 


Background 


9.36 


7.83 


8.59 


10.82 


Polarization 


6.11 


4.72 


4.67 


7.19 


Mc Statistics 


3.45 


1.68 


2.72 


2.27 


Total 


18.42 


15.95 


16.13 


19.14 



Table 4: Breakdown of contributions to the systematic errors for the Xci K* channels. All values 
are expressed relative to the measured branching fractions, in percent. 







K+ -K- 




K+ ^0 


B counting 


1.1 


1.1 


1.1 


1.1 


Tracking 


2.6 


5.2 


3.9 


3.9 




1.1 




1.1 




Neutral 


9.4 


2.5 


2.8 


8.3 


Second BF 


10.4 


10.4 


10.4 


10.4 


PID 


0.35 


1.64 


0.85 


1.23 


Background 


30.08 


16.76 


24.25 


246.53 


Polarization 


8.27 


5.86 


6.81 


8.61 


MC statistics 


2.06 


1.38 


1.61 


1.80 


Total 


34.40 


21.50 


27.78 


247.08 



channel-dependent correction. For the particle identification efficiency correction, we have assigned 
a systematic error equal to half of the correction. The overall selection efficiency depends on 
the angular distribution used in the simulation for the decay. It can be written as e = a + | Aop6, 
where a and h are functions of the K* helicity angle, a = 3/4/ (1 — cos^ Sk*)^{Ok*) sm{9K*)d6K* and 
h= 3/4 J(3 cos^ ^if* — l)e(0i^*) sin(0/^»)(i0;^* , and |^ol is the (unknown) fraction of the longitudinal 
K* polarization |^. We estimate the uncertainty due to our ignorance on the value of |j4o| and 
derive an associated systematical error varying from 3.4 to 8.6%, depending on the channel. The 
systematic error due to the finite size of the Monte Carlo statistics varies from 1.23 to 3.45%. 
Interference effects between K* events described by a P-wave and non-resonant events described 
by an S-wave have been considered ^U]. The interference term is proportional to the fraction of 
non-resonant events with respect to the number of signal events ^01 ■ However, this fraction is small 
for all channels. Furthermore, a large systematic uncertainty (see below) has been assigned to the 
number of non-resonant events. Thus, no additional systematic uncertainty due to interference 
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Table 5: Breakdown of contributions to the systematic errors for the K'^ and channels. All 
values are expressed relative to the measured branching fractions, in percent. 





e+e M^A* 


J/i) K+ 
e+e~ /Lt^M^ 


e+e /^^/^ 


^j{2S) K+ 
e+e^ jJL'^ 


e+e A^^M 


Xci K+ 

e+e fJ-^ l-J- 


B counting 
Tracking 

PID 

Second BF 
Background 
MC statistics 


1.1 
2.6 
0.7 

0.20 1.55 
1.70 1.71 
0.22 0.08 
0.56 0.54 


1.1 
3.9 

1.04 1.55 
1.69 1.70 
0.17 0.03 
1.07 1.01 


1.1 
2.6 
1.0 

0.22 2.41 
4.11 10.96 
0.41 0.27 
1.36 1.31 


1.1 
3.9 

0.74 2.08 
4.11 10.96 
0.57 0.17 
1.88 1.81 


1.1 
2.6 
0.9 

0.74 1.94 
10.27 10.26 
6.72 2.96 
1.14 1.11 


1.1 
3.9 

0.82 2.67 
10.27 10.27 
1.97 1.84 
1.54 1.50 


Total 


3.46 3.75 


4.64 4.77 


5.28 11.69 


6.14 12.01 


12.70 11.30 


11.40 11.665 



effects has been introduced. 

In the default fit for the determination of the combinatorial background, the shape parameter 
of the Argus function is not constrained. To determine a systematic error, a second fit with 
the shape parameter of the Argus function fixed to the value obtained from fitting the data in 
the AE sideband region was performed. We have taken as the systematic uncertainty on the 
combinatorial background 50% of the difference between the combinatorial background contribution 
obtained from the default fit and from the second fit. For the feed-across component to the 
peaking background we have assigned as the systematic error, the uncertainty of the corresponding 
branching fractions, taken from _8^. For the inclusive charmonium contribution to the peaking 
background, we have assigned a 50% error, accounting for the poor knowledge of the branching 
fractions of the contributing decay modes. 

The systematic uncertainties for all modes are listed in Tables El El El and El 

5 Physics results 

The branching fractions that have been measured separately for the J/ip — > e~^e~ and J/tp 
IJL^ IJi~ decay modes were found to be in good agreement. They have therefore been combined. 
Furthermore, for K* channels, the branching fractions from the two neutral sub-modes tt^ and 
i^+ 7r~ have been averaged together, and the branching fractions from the two charged sub-modes 
7r+ and i^+ vr*^ have been averaged together as well. The branching fraction measurements are 
summarized in Table El 

From these measurements, we have determined the ratios of charged to neutral branching 
fractions. We have assumed a value of one for the charged to neutral B meson production rate 
at the T(45'). The results are presented in Table [71 The systematic uncertainties of the ratios 
have been determined by taking into account the correlations of the errors between the branching 
fractions. 

Combining all the measurements, we obtain: 



B(B~^ charmonium K(*)+) . , 

^ = 1.14 ±0.02 ±0.03 (2 

o(B^ — >■ charmonium K"-*^") 

15 



Table 6: Measured branching fractions for exclusive decays of B mesons to charmonium and kaon 
final states. The first error is statistical and the second systematic. 



Channel 


Branching fraction (xlO 


- 


J/ij) K*^ 


12.92±0.25±0.75 


B+ - 




14.34±0.36±0.94 


B+ 


J/ih K+ 

IT 


10.55±0.15±0.48 


B^ - 


J/ilj K° 


8.73±0.23±0.30 


B^ - 


i^{2S) K*^ 


6.65±0.57±1.00 


B+ 


V(25) K*+ 


6.03±0.85±0.91 


B+ 


V(25) K+ 


6.31±0.33±0.44 


i^o - 


V(25) 


6.60±0.60±0.46 


B^ - 


Xci K*o 


3.19±0.37±0.64 


B+ 


- Xci K*+ 


2.89±0.69±0.93 


B ^ 


► Xci K+ 


5.72±0.24±0.64 


B^ - 


Xci 


4.56±0.39±0.51 



Table 7: Results for ratios of charged to neutral braching fractios. The first error is statistical and 
the second systematic. 



Ratio 






Result 


B{B+ - 


J/i)K+)/B{B'' ^ 




1.21±0.04±0.04 


B{B+ - 


'tp{2S)K+)/B{B° 


yj{2S)K^) 


0.95±0.10±0.03 


B{B+ - 


XciK+)/B{B^ ^ 


XciK^) 


1.25±0.12±0.07 


B{B+ - 


J/i;K*+)/B{B^ - 


J/^K*^) 


1.11±0.04±0.08 


B{B+ - 


i;{2S)K*+)/B{B^ ^{2S)K*^) 


0.91±0.15±0.11 


B{B+ - 


XciK*+)/B{B^ - 


> XciK*^) 


0.91±0.24±0.31 



Assuming isospin invariance in the B — > charmonium K {K*) decays we can compute our own 
value for the charged to neutral B meson production. Using the ratio of the charged to neutral B 
meson lifetimes t^+It^o = 1.086 ± 0.017 ^5), we obtain: 



r(r(4S) ^ MO) ^ ' 

We also determine the ratio of branching fractions for a vector to a pseudo-scalar light meson: 
B{B^ ipK*°)/B{B° i^K^) and B{B+ i;K*+)/B{B+ tpK+) for the three charmonium 
states = J/ip , ip{2S) and Xci- The results are presented in Table |H1 For each of the charmonium 
states, we also present the average of the charged and neutral measurements. 

Finally, charge asymmetries have been measured. The branching fractions for positively and 
negatively charged B mesons have been determined using the method described above. The selec- 
tion efficiencies have been determined separately. 



16 



Table 8: Results for ratio of the branching fractions for a vector (K*) versus pseudoscalar (K) 
light meson. The first error is statistical and the second systematic. 



Ratio Result 

B{B^ J/^K*'^)/B{B'^ J/^jK"^) 1.48±0.05±0.07 

bIb+ Jj^ K*+)/B{B+ J/^ K+) 1.36±0.04±0.08 

B{B J/%l)K*)/B{B ^ J/i)K) 1.42±0.03±0.05 

^ i){2S)K*'^)/B{B'' ij{2S)K'') 1.01±0.13±0.09 

B{B+ i^{2S)K*+)/B{B+ V(25)K+) 0.96±0.14±0.09 

B{B %b{2S)K*)/B{B ^{2S)K) 0.99±0.10±0.06 

BiB"^ XciK*^)/B{B" XciK") 0.70±0.10±0.12 

B{B+ XciK*+)/B{B+ XciK+) 0.51±0.12±0.15 

B{B XciK*)/b[b XciK) 0.62±0.08±0.09 



B{B+ 


J/^K+) 


- B{B- - 


^ J/i)K-) 


B{B+ 


J/^K+) 


+ B{B~ - 


■> J/^l}K~) 


B{B+ 


J/tpK*+ 


)-B{B- 




B{B+ 


J/tpK*+ 


)+B{B~ 


J/^K*-) 


B{B+ - 


ij{2S)K+ 


)-B{B' 


-^i,{2S)K-) 


B{B+ - 


^{2S)K+ 


)+B{B- 


^i^{2S)K-) 


B{B+ 


ij{2S)K*+) 


- B{B- - 


■> i){2S)K*-) 


B{B+ 


ipi2S)K*+) 


+ B{B~ - 


■> i){2S)K*-) 


B{B 


^ ^ XciK+ 


) - B{B- 


- XciK~) 


B{B 


^ ^ XciK+) + B{B- 


- XciK~) 


B{B+ 




- B{B- - 


- XciK*-) 


B{B+ 




+ B{B- - 


- XciK*-) 



0.029 ±0.014 ±0.010 



0.045 ±0.025 ±0.011 



0.059 ± 0.051 ± 0.021 



-0.063 ±0.137 ±0.050 



0.011 ±0.042 ±0.017 



-0.403 ± 0.309 ± 0.237 



6 Summary 

We have presented preliminary results of branching fraction measurements of exclusive B decays 
to charmonium and K or K* . The charmonium mesons considered were J/ip , '4>{2S) and Xci- Our 
results for J/ip and ip{2S) are in good agreement with previous measurements |S] with comparable 
or superior precision. Our Xci results have much better precision and the B~^ Xci K*^ branching 
fraction was measured for the first time. Assuming isopin invariance, we find the ratio of charged 
to neutral B meson production on the T(4S) resonance to be compatible with unity within two 
standard deviations. No direct CP violation has been observed from the measurements of charge 
asymmetries as we found them to be compatible with zero. 
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